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Abstract 

The light cone QCD sum rules method is used to calculate the transition form 
factor for the radiative (j) — > /07 decay, assuming that the quark content of the /o 
meson is pure ss state. The branching ratio is estimated to be B{4) — > /07) = 
3.5 X (1 lb 0.3) X 10~^. A comparison of our prediction on branching ratio with the 
theoretical results and experimental data existing in literature is presented. 
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1 Introduction 



According to the quark model, mesons are interpreted as pure qq states. Scalar mesons 
constitute a remarkable exception to this systematization and their nature is not well es- 
tablished yet [l]-[4]. 

In the naive qq picture, one can treat the isoscalar /o(980) either as the meson that 
exists mostly as nonstrange and almost degenerate with the isovector ao(980) or as mainly 
ss, in analogy to the pure ss vector meson 0(1020). 

In order to understand the content of the fo meson several alternatives have been sug- 
gested, such as, the analysis of the fo — * 27 decay [5, 6]; study of the ratio r(ao — > 
/o7)/r(0 /07) [7]. However, among these, the (0 — > /07) decay occupies a special place, 
since the branching ratio expected of this decay, is essentially dependent on the content of 
fo. For example i3(0 /07) is as high as ~ 10^^ if it were composed of qqqq, and ~ 10~^ 
if /o were a pure ss state. 

It has been known for a long time that /o(980) couples significantly through its ss 
content, from its detection as a peak in the J/ip 0/o [8] and Ds nfo [9] decays, as 
discussed in [10] and [11] (see also [12]). For this reason, in this work we assume that quark 
content of both (p and fo mesons are pure ss. In the present paper we analyze the radiative 
~^ /o7 decay in framework of the hght cone QCD sum rules (about hght cone QCD 
sum rules and its applications, see for example [13]). Note also that the — >■ /07 decay is 
analyzed in framework of the 3-point sum rules in [14]. In order to calculate the transition 
form factor describing the (j) — > /07 decay in light cone QCD sum rules, we consider the 
following correlator 

n^^ij d^xe"^^ {o\T{r{x)j^{o)}\o)^ , (1) 

where J* = ss and = s'-f^s are interpolating currents for fo and mesons, respectively, 
and 7 is the background electromagnetic field (for more about external field technique in 
QCD see [15, 16]). 

The physical part of the correlator can be obtained by inserting a complete set of one 
meson states into the correlator. 



^ _ ^ {0\J%x)\fo{p)) ifoipMpi)), (0(pi)| J^^(0)|0) 
[p^ - m%) [pi - m; J 

where and /o are the quantum numbers and pi — p+q with q being the photon momentum. 
The matrix element ^0(pi) -'^(O) 0^ in Eq. (1) is defined as 

(0(pi)|^^(O)|o) = m^/<^£^, (3) 

where e:^ is the meson polarization vector. The coupling of the /o(980) to the scalar 
current J* = ss is defined in terms of a constant Xf 

{Q\J'\fo{p))=mf,Xf . (4) 
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The relevant matrix element describing the transition (j) ^ fo induced by an external 
electromagnetic current can be parametrized in the following form: 



(5) 



where e is the photon polarization and we have used {sq) = 0. From gauge invariance we 
have 



(6) 



and since the photon is real in the decay under consideration, we need the values of the 
form factors only at the point q'^ — 0. Using Eq. (6) the matrix element (/o|0)^ takes the 
following gauge invariant form, 



{fo\ci>)^^ee^F,{0) {p,q)sf, - {e'^q)p 



Using Eqs. (l)-(4) and (7), for the phenomenological part of the correlator we have 



rf- = eFi(OK -{Piq)9f..+Piuq^ 



(f - m%) (p\ - mfj 



(7) 
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In order to construct the sum rule, calculation of the correlator from QCD side (theoretical 
part) is needed. Prom Eq. (1) we get 







(9) 



where Sg is the full propagator of the strange quark (see below). Theoretical part of 
the correlator contains two pieces, perturbative and nonperturbative. Perturbative part 
corresponds to the case when photon is radiated from the freely propagating quarks. Its 
expression can be obtained by making the following replacement in each one of the quark 
propagators in Eq. (9) 



Q oh 



ah 

al3 ' 



(10) 



where the Fock-Schwingcr gauge x'^A^{x) = is used and S'/''^'^ is the free s-quark propa- 
gator S{^'^'^{x) = i ^1 (^11^ x^) and the remaining one is the full quark propagator. 

The nonperturbative piece of the theoretical part can be obtained from Eq. (9) by 
replacing each one of the propagators with 



(11) 



where is the full set of Dirac matrices and sum over is implied and the other quark 
propagator is the full propagator, involving pertubative and nonperturbative contributions. 
In order to calculate perturbative and nonperturbative parts to the correlator function (1), 

expression of the s-quark propagator in external field is needed. 

The complete light cone expansion of the light quark propagator in external field is 
presented in [16]. The propagator receives contributions from the nonlocal operators qGq^ 
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qGGq^ qqqq, where G is the gluon field strength tensor. In the present work we consider 
operators with only one gluon field and neglect terms with two gluons qGGq, and four 
quarks qqqq and formal neglect of these these terms can be justified on the basis of an 
expansion in conformal spin [17]. In this approximation full propagator of the s-quark is 
given as 

i ^ (ss) 



Ss{x) 



X 



2tt'^x^ 



12 



,^+16'"°, 



ims(ss) , irninis o 
H — — ^ X 



48 



2^32 



dv 



:G,. 



:VX''G,, 



(12) 



It follows from Eqs. (11) and (9) that in calculating the QCD part of the correlator, as 
is generally the case, we arc left with the matrix elements of the gauge invariant nonlocal 
operators, sandwiched in between the photon and the vacuum states (7(g) |sAjs| 0). These 
matrix elements define the light cone photon wave functions. The photon wave functions 
up to twist-4 are [17, 18] 



(7(g) \q{x)a^^q{0)\0) = ieeg{qq) ( due"^''[ {e^q^ - £^g^) [x0(w) + x'^(^gi{u) - g2{u)) 

J 



+ 



(gx)(e^x^ 



+ 



{ex){Xf,q^ - x^qf,)jg2{u)} , 



(13) 



(7(g) |g(a;)7^75g 



0) = —egea(3p<7£^q''x" 



(14) 



The path-ordered gauge factor Vexp(igs Jq dux^A^{ux)^ is emitted since the Schwinger- 
Fock gauge a;'^A^(x) = is used. The functions (f){u), ip{u) are the leading twist-2 photon 
wave functions, while gi{u) and g2{u) are the twist-4 photon wave functions. Note that 
twist -3 photon wave functions are neglected in the calculations, since their contributions are 
small and change the result by 5%. In Eq. (13) x is the magnetic susceptibility of the quark 
condensate and Sq is the quark charge. The theoretical part is obtained by substituting 
photon wave functions and expression for the s-quark propagators into Eq. (9). The sum 
rules is obtained by equating the phenomenological and theoretical parts of the correlator. 
In order to suppress higher states and continuum contribution (for more details see [19, 20]) 
double Borel transformations of the variables p\ — and P2 — {p + q)^ are performed on 
both sides of the correlator, after which the following sum rule is obtained 



Fi(0) 



2x{ss)(p{uo) - 



2^ 



M^Eo{sl/M^ 



+ ^(^^) - l92gi{uo) + ms4>{uo){ss) 



+ 



3ms 
2^ 



+ In^ j Eo{so/M') + M'f{so/M' 



(15) 



where Sq is the continuum threshold 

Eo{so/M^) = 1 - e-^°/^' 



Jo 
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which have been used to subtract continuum, and 



Ml + M| 



-2 ' 



Ml + M| 



■2 ' 



where and M| are the Borel parameters in and /o channels, respectively, A is the 
QCD scale parameter and 7e is the Euler constant. Since the masses of (f) and /o are very 
close to each other we will set Ml — M| = 2M^, obviously from which it follows that 



It is clear from Eq. (14) that the values of A/g and are needed in order to determine 
F(0). The coupling of the /o(980) to the scalar ss current is determined by the constant 
A/o and in the two-point QCD sum rules its value is found to be A/g = (0.18 ±0.0015) GeV 
[14]. In further numerical analysis we will use = 0.234 GeV which is obtained from the 
experimental analysis of the — > e^e" decay [21]. 

Having the values of A/,, and our next and final attempt is the calculation of transition 
form factor i^i(O). As we can easily see from Eq. (15) the main input parameters of the 
light cone QCD sum rules is the photon wave function. It is known that the leading photon 
wave function receive only small corrections from the higher conformal spin [17, 19, 22], so 
that they do not deviate much from the asymptotic form. The photon wave functions we 
use in our numerical analysis are given as 



Furthermore, the values of the input parameters that are used in the numerical calculations 
are: / = 0.028 GeV"^ , x = -4.4 GeV'^ [23] (in [24] this quantity is predicted to have 
the value x = -3.3 GeV-^), {ss(l GeV)) = -0.8 x (0.243)''^ GeV^ and the QCD scale 
parameter is taken as A = 0.2 GeV. The strange quark mass is chosen in the range 
nis — 0.125 0.16 GeV, obtained in the QCD sum rules approach [25]. The masses of 
the (f) and /o mesons are — 1.02 GeV, mf^ — 0.98 GeV. The transition form factor 
is a physical quantity and therefore it must be independent of the auxiliary continuum 
threshold Sq and and the Borel mass M^ parameters. So our main concern is to find a 
region where the transition form factor -Fi(O) is practically independent of the parameters 
So and M^. For this aim in Fig. (1) we present the dependence of the transition form factor 
Fi(0) on the Borel parameter at three different values of the continuum threshold: 
So = 2.0 GeV^, 2.2 GeV^ and 2.4 GeV^. It follows from this figure that for the choice 
of the continuum thresholds in the above-mentioned range, the variation of the result on 
the transition form factor Fi(0) is about 10%. In other words, we can conclude that Fi(0) 
is practically independent of the continuum threshold. Furthermore we observe that when 
1.4 < M^ < 2.0 GeV'^, -Fi(O) is quite stable with respect to the variations of the Borel 
parameter M^. As a result, one can directly read from this figure 



Uo = 1/2. 



(f){u) = 6u{l — u) , 
iP{u) = 1 , 




Fi(0) = (3.25 ±0.20) GeV 



-1 
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where the resulting error is due to the variations in sq and M^. The other sources of 
errors contributing to the numerical analysis of the transition form factor come from the 
strange quark mass and the uncertainties in values of various condensates. Hence, our final 
prediction on the transition form factor is 

Fi(0) = (3.25 ± 0.50) GeV-^ . (16) 

Using the matrix element (7) for the decay width of the considered process, we obtain 

- hi) = a \FMf ^ \^J^'^ ■ (17) 

Using the experimental value rtot(0) = 4.458 MeV [21], and Eqs. (16) and (17), we get for 
the branching ratio 

-6(0 ^ /o7) = 3.5 X (1.0 ± 0.3) X 10"^ . (18) 

Our result on the branching ratio is obtained under the assumption that /o meson is 
represented as a pure ss component. How does the result change if we assume that and 
/o mesons can be represented as a mixing of ss and nn = {uu + dd)/\/2 state, i.e., 

(f) — cos ass -\- sin a fin , 
/o = sin (5ss-\- cos (5nnl 

Analysis of the process ti^'~^ and combined analysis of the /o7 and /o — > 27 decays 
show that \a\ < 4° and two solutions are found for /3, i.e., P = -48° ± 6° or /3 = 85° ± 5°, 
respectively (see for example [5]). In other words, quark content of (p meson is pure ss 
state, while in /o meson there might be sizable nn component. Obviously, when -Fi(O) is 
calculated from QCD side, only sin f3 ss component operates (see Eq. (1)) and therefore the 
decay width r(0 — > /07), and hence the corresponding branching ratio, contains an extra 
factor sin^/3. li P — 85° ± 5°, then prediction for the branching ratio given in Eq. (18) is 
practically unchanged, but when /3 = —48° ± 6° B(0 — > /07) decreases by about a factor of 
2. 

Finally let us compare our prediction on branching ratio with the existing theoretical 
results and experimental data in the literature. Obviously, our result is slightly larger 
compare to the 3-point QCD sum rule result which predicts B{(p fol) — (2.7±1.1) x 10"^ 
[14] , and approximately three times larger compared to the prediction of the spectral QCD 
sum rules and chiral unitary approaches, whose predictions are S(0 fo'~f) = 1.3 x 10~^ 
[26] and B{(f) /07) = 1.6 x 10^^ [27], respectively. It is interesting to note that this value 
of the branching ratio is closer to our prediction when the mixing angle is chosen to be 
P — —48° ±6°. Our result, which is given in Eq. (18), is larger compared to the predictions 
of [28, 7], whose results are B(0 ^ /07) = 1.9 x 10"^ [28], and B(0 ^ /07) = 1.35 x 10"^ 
[7], respectively. 

As the final words we would like to point out that our prediction given in Eq. (18), is in a 
very good agreement with the existing experimental result i3(0 — >■ /07) — (3.4 ± 1.1) x 10~"^ 
[21]. 
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Figure 1: The dependence of the transition form factor Fi{0) for the radiative — > /o7 
decay on at three different values of the continuum threshold Sq — 2.0 GeV"^, 2.2 GeV"^ 
and 2.4 GeV^. 
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